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We give a pedagogical review of implications of chiral symmetry in QCD. First, 
we briefly discuss classical textbook subjects such as the axial anomaly, sponta- 
neous breaking of the flavor-nonsinglet chiral symmetry, formation of light pseudo- 
Goldstone particles, and their effective interactions. Then we proceed to other is- 
sues. We explain in some detail a recent discovery how to circumvent the Nielsen- 
Ninomiya's theorem and implement chirally symmetric fermions on the lattice. 
We touch upon such classical issues as the Vafa-Witten's theorem and 't Hooft's 
anomaly matching conditions. We derive a set of exact theorems concerning the 
dynamics of the theory in a finite Euclidean volume and the behavior of the Dirac 
spectral density. Finally, we discuss an imaginary world with a nonzero value of 
the vacuum angle 9. 
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1 Classic Sagas 

1.1 Basics 

When a field theory enjoys a symmetry, it is always a fortunate circumstance 
because the presence of the symmetry always simplifies analysis and, more 
often than not, allows one to obtain exact or semi-exact results. The richer 
is the symmetry, the stronger are the results thus obtained. Sometimes, the 
presence of a rich enough symmetry allows one to solve the theory (the best 
example is the exact results for the spectrum .of the M = 2 supersymmetric 
Yang-Mills theory due to Seiberg and WittenEl). The symmetry of QCD is 
not so rich and the theory is not exactly solvable. Still, the presence of light 
quarks such that the theory enjoys chiral symmetry in the chiral limit, when the 
quarks become massless, allows one to extract a lot of consequences concerning 
QCD dynamics, with some of them having the status of exact theorems. In 
this mini-review we will concentrate on some less known and/or comparatively 
recent results obtained in this direction. To make the discussion coherent, we 
will review some well-known facts as well. 
The Lagrangian of QCD reads 

Cqcd = -^T?{F, v Fn + J2 hdf- ™Mf + 4^ F ^ a . « 
9 f=i 

where = ^e^vapFap, the physical value of the parameter 9 is very small 
9 < 1CP 9 , and rrif are the quark masses. The u, d, and s quarks are relatively 
light, with masses m u w 4 MeV, w 7 MeV, and m s w 150 MeV. m U)( j are 
especially small. It makes sense to consider the chiral limit when Nf = 2 or 
Nj — 3 quarks become massless. In this limit, the Lagrangian ([l]) is invariant 
under the transformations 



8^ f =ia A [t A i>]f (2) 



and 



<% = */3a7 5 [^]/ (3) 

where t A (A = 0, 1, . . . , Nj — 1) are the generators of the flavor XJ(Nf) group. 
The symmetry (||) is the ordinary isotopic symmetry. It is still present even 
if the quarks are endowed with a mass (of the same magnitude for all fla- 
vors) . The symmetry (^|) holds only in the massless theory. The corresponding 
Noether currents are 

if) A = ^v^, (r 5 ) A = v^VtV- (4) 

They are conserved upon applying the classical equations of motion. 
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1.2 Singlet Axial Anomaly 

Let us discuss first the singlet axial symmetry (with t A = 1). An important 
fact is that this symmetry exists only in the classical case. The full quantum 
path integral is not invariant under the transformations 

Sip = iaj 5 ^, StJ) — iofip^ 5 , (5) 

where the flavor index is now omitted. This explicit symmetry breaking due to 
quantum effects can be presented as an operator identity involving an anoma- 
lous divergence, 

d,f 5 = -^e a ^Tr{F a0 F^} . (6) 

There are many ways to derive and understand this relation. Historically, it 
was first derived by purely diagrammatic methods. The amplitude 

S a % uX (k,q) = /d 4 xdV fcx e^(T{j-(x)^(0)jt(y)})o, (7) 



(ja — $f 7/iV") where a, b are color indices) cannot be rendered transverse 
with respect to both vector and axial indices such that k^T^ v \ — and 
<l\Tfi U x = due to ultraviolet divergences in the anomalous triangle graph. 
We will comment more on this later when we will be discussing 't Hooft's 
self-consistency conditions. 

The anomaly relation (^|) can also be derived as an operator identity. To 
this end, one should carefully define the operator of the axial current as 

f 5 = Vm^(x + e)^expUj + A a (y)dy a ^{x) (8) 

and perform the limit, averaging over the directions of and making use of 
the operator equations of motion (see e.g. Ref. 2 for more details.). 

Let us dwell here on the third and the most elegant way to derive the 
anomaly relation (|^) which uses path integral methods. The anomaly appears 
due to the necessity to regularize the theory in the ultraviolet. The most 
politically correct approach would be to study a path integral regularized by 
a lattice, and we will do so in Sec. 2. Let us, first, describe a more standard 
and habitual finite mode regularization. 

Consider an Euclidean path integral for the partition function in QCD 
with one massless quark flavor. The fermionic part of the integral is 



Yl dip(x)dt/j{x) expli / d i xi/jp E ip 
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which formally coincides with the determinant of the Euclidean Dirac operator 
—Hf> E — —ij^V^, where 7^ are Euclidean 7 matrices. The latter are anti- 
Hermitean and satisfy 



From now on, we will stay in the Euclidean space and will suppress the super- 
script "E" ioxTp and 7^. 

Let us assume that the theory is somehow regularized in the infrared so 
that the spectrum of the operator^ is discrete. Note now that the spectrum of 
the massless Euclidean Dirac operator enjoys the following symmetry: for any 
eigenfunction Uk of the operator Tp with an eigenvalue , the function u' k — 
j 5 Uk is also an eigenfunction with the eigenvalue —Xk tjp and 7 s anticommute) . 
Therefore, 

1 

det||-^ + m|| = JJ(-iAfc+m) = m q ]J{\ 2 k + m 2 ) , (9) 

k k 

where q is the number of possible exact zero modes and the product Y\' k runs 
over the nonvanishing eugenvalues. We see that the expression ([)]) is real, and 
the Euclidean partition function is real, too. 
Let us expand ip{x), ip(x) as 

{ - t (10) 

{ i>{x) = J2k c kul{x) 

where {uk{x)} is a a set of eigenfunctions ofTp forming a complete basis in the 
corresponding Hilbert space. Then 

Y[dtp(x)dip(x) = Y\dc k dc k - (11) 

x k 

Suppose the field variables are transformed by an infinitesimal global chiral 
transformation (|f])|] Now, ip' — ip + Sip and $ — "tp + S'tp can be again 
expanded in the series (|l0|). The new expansion coefficients are related to the 
old ones, 

m ' m 

c' k = c k + ia^2c m J d 4 x u^x)^ u k {x) = ^ c m (S km + iaA mk ) ■ (12) 



a We keep the Euclidean ip and ip independent and are not bothered by the fact that Sip ^ 

W) f . 



The point is that the transformation ( |l^ ) has a nonzero Jacobian. We have 
J|dcfe'd4 = J~ 2 Y[ dc k dc k , (13) 



where 



J = dct(l + iaA) s» exp 



(14) 



Or, in other words, 

In J = ia I AVti|(i)7 5 !i);(i) + o(a) . (15) 



k 



One's first (wrong!) impression might be that J d 4 x^2 k u k (x)j 5 Uk(x) is just 
zero. Indeed, as was already mentioned above, the function u' k = j 5 Uk is also 
an eigenfunction of the Dirac operator with the eigenvalue —Xk- If Afc 7^ 0, 
Uk(x) and 7 5 u/c(x) thereby represent different eigenfunctions and the integral 
J d 4 x u\.{x)~f b Uk{x) vanishes. 

A nonzero value of ( |l5| ) is due to the fact that, for intricate enough topo- 
logically nontrivial gauge fields, the spectrum of the Dirac operator involves 
some number of exact zero modes, for which 7 5 uo(x) = =F 1/0(2;) (depending on 
whether the modes are left-handed or right-handed) , and their contribution to 



the sum (15) is responsible for the whole effect. A famous theorem of Atiyah 
and Singer dictates 

J d 4 xJ2u[(xh 5 Mx) = r$ - n f = 1 > (16) 

k 

where nf^ R is the number of the left-handed (right-handed) zero modes and 
q is the topological charge of the gauge field configuration, 

q = ik* / d4xTr i F ^ p ^}- 

Substituting (|l^) in Eqs. (|ll| |l3|), we now see that the change of the 
measure under the chiral transformation can be presented as a shift of the 
effective action 

SS = Y^2 e «/3^ / d 4 xTr{F a/3 F^ .} (18) 
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In other words, a global singlet chiral transfomation is equivalent to leaving 
the fermionic fields intact, but shifting instead the parameter 6 in the original 
theory §), 

6^9 + 2a. 

1.3 Spontaneous Breaking of Nonsinglet Chiral Symmetry 

Consider now the whole set of symmetries @ , ([}]) . It is convenient to introduce 

4l,r= ^(1T1 5 )^, rpL,R= ^(1±7 5 ) (19) 

and rewrite (0), (||) as 

ipL -» V L i> L , i> R -» V R ^ R , (20) 

where Vl and Vr are two different U(iV/) matrices. The singlet axial transfor- 
mations with Vl = Vr = e 1 ^ are anomalous by the same token as in the theory 
with a single quark flavor. Therefore, the true fermionic symmetry group of 
massless QCD is 

Q = SU L (JV» x SU fl (JV» x U v (l) ■ (21) 

A fundamental experimental fact is that the symmetry ( pl| ) is actually sponta- 
neously broken, which means that the vacuum state is not invariant under the 
action of the group Q. The symmetry Q is, however, not broken completely. 
The vacuum is still invariant under transformations with Vl = Vr, generated 
by the vector isotopic current.^] 

Thus, the pattern of breaking is 

SUi(iV» x SU H (JV» -► SVv(N f ) . (22) 

The nonvanishing vacuum expectation values 

E/» = (ij/J R ) (23) 

are the order parameters of the spontaneously broken axial symmetry. The 
matrix is referred to as the quark condensate matrix. 

Non-breaking of the vector symmetry implies that the matrix order pa- 
rameter ( p3| ) can be cast in the form 

E /9 = ~E5 /9 (24) 



'An exact theoremQ that the vector symmetry cannot be broken spontaneously in QCD will 
be proven in Sec. 3. 
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by group transformations (|20|). This means that the general condensate matrix 
is a unitary SU (Nf) matrix multiplied by E. 

In general, £ could be any complex number. It can be made real by a global 
Uyi(l) rotation which, according to Eq. ( fl8|) (with the factor Nf), amounts to 
a shift of the vacuum angle 9. For massless quarks, physics does not depend 
on the phase of £ and, thereby, on 9. It is convenient then to choose £ real 
and positive and 9 = 0. From experiment, we know that E ss (240 MeV) 3 /2 
with about 30% uncertainty (this value refers to some particular normalization 
point [i ~ 0.5 GeV on which the operator ifiip and its vacuum expectation value 
depend) . 

By Goldstone's theorem, spontaneous breaking of a global continuous sym- 
metry leads to the appearance of purely massless Goldstone bosons. Their 
number coincides with the number of broken generators, which is Nj — 1 in 
our case. As it is the axial symmetry which is broken, the Goldstone parti- 
cles are pseudoscalars. They are nothing but pions for Nf = 2 or the octet 
(<ir,K,K,ri) for Nf = 3. 

It is a fundamental and important fact that spontaneous breaking of con- 
tinuous symmetries not only creates massless Goldstone particles, but also fixes 
the interactions of the latter at low energies. To see this, let us first recall that 
the Goldstone field describes fluctuations of the order parameter, 

E fg -> Z fg {x) = -EU(x) 

with U(x) € SU(Nf). It is convenient to express U(x) as an exponential 

, s (2i6 a (x)t a } , s 

U(x) = exp ^ > , (25) 



Fir 

where 4> a (x) are the physical meson fields [so that (f> a (x) = corresponds to 
the vacuum (p4|)] and F n is a constant of dimension of mass. 

The Goldstone particles are massless whereas all other states in the phys- 
ical spectrum have nonzero mass. Therefore, we are in the Born-Oppenheimer 
situation: there are two distinct energy scales and one can write down an effec- 
tive Lagrangian depending only on slow Goldstone fields with the fast degrees 
of freedom corresponding to all other particles being integrated out. 

To fix the exact form of this Lagrangian, note that the transformations 
( f2p| ) are realized at the level of the effective Lagrangian as U — * VzUVl. Any 
scalar function depending on U and invariant under this symmetry is also a 
function of U^U = 1, i.e. it is just a constant. There is only one invariant 
structure involving two derivatives, namely 

4ff = ^fTv{d^Ud^} . (26) 
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Take Nf = 2. The perturbative expansion of Eq.(p6|) in powers of 4> reads 

4ff = l(d,r) 2 + ^ [ow*) 2 - (rr)(d^ b ) 2 ] + ... (27) 

Also for Nf > 3, we have, on top of the standard kinetic term, a quartic term 
involving two derivatives, with somewhat more complicated group structure. 
We see that the symmetry dictates rather specific interactions between the pi- 
ons. They do not interact at the s-wave level, which means that all amplitudes 
vanish at zero momenta, but the strength of interaction grows rapidly with 
energy. 

Equation (|2(]) describes the effective chiral Lagrangian of massless Gold- 
stone bosons at leading order. The first corrections involve 4 derivatives and 
there arc 3 different invariant functions 

4ff = iiTrj^C/d^} 2 +L 2 Tr{d f ,Ud„U t }Tr{d f *Ud"U t } 

+ L 3 Tr{d ll Udi J -U^d v Ud l/ U^} (28) 

(only 2 linearly independent structures are left for Nj = 2). 

The relevant Born-Oppenheimer expansion parameter is K c hir ~ p char JF V . 
When K c hir ~ 1 (in practice, one should rather take K c hir ~ 27r), the Born- 
Oppenheimer approximation as well as the whole effective Lagrangian ap- 
proach breaks down, and non-Goldstone degrees of freedom become impor- 
tant. The physical meaning of F w is thus clarified. It characterizes the gap in 
the spectrum and sets a scale below which massive degrees of freedom can be 
disregarded. 

Let us discuss the actual world now. The Lagrangian of real QCD (@) is 
not invariant under the axial symmetry transformations just because quarks 
have nonzero masses. The symmetry ( |2l| ) is still very much relevant to QCD 
because some of the quarks happen to be very light. 

For Nf = 2, spontaneous breaking of an exact SUl(2) x SUi{(2) symmetry 
would lead to the existence of 3 strictly massless pions. As the symmetry is 
not quite exact, the pions have a small mass. However, their mass M goes to 
zero in the chiral limit m u ^ — > 0. Indeed, trading the mass term 

- m u uu - m d dd = m u (u L u R + u R u L ) + m d {d L d R + d R d L ) (29) 

in the QCD Lagrangian for the contribution 

C[p = S Re [Tr{MU*}] (30) 
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(M is the quark mass matrix which is chosen here in the form M = 
diag(m u , m<j) with real m u ,md) in the effective chiral Lagrangian^ and ex- 
panding ( |30| ) in 4> a , we obtain the G ell- Mann- Oakes-Renner relation 

FlM 2 = (m u + m d )E + 0(m 2 q ) . (31) 

The constant F„ appears also in the matrix element 

(vac|A+|vr) p =iV2F„pl 

of the axial current — dj^^u and determines the charged pion decay rate. 
Experimentally, F w « 93 MeV. 

The effective chiral Lagrangian involving the lowest ( p6f ) , (|3^) , and higher 
order terms (in addition to the terms of higher order in derivatives as in 
Eq. (p£j|), there are also terms of higher order in mass) can be and was used to 
calculate the amplitudes of different processes involving pseudoscalar mesons at 
low energies. The corresponding technique is called chiral perturbation theorya 

2 Quarks on Euclidean Lattice 

The material of the previous section is well known and can be found in many 
textbooks. We are in a position now to discuss less known issues and will start 
with the question of how to put quarks on the lattice. This issue looked very 
confusing for a long time, and has been clarified only recenly. 

2.1 The Nielsen-Ninomiya's No-go Theorem 

The only way to define what quantum field theory in nonperturbative regime 
really means is to introduce a lattice ultraviolet regularization so that the 
symbol of path integral is defined as a limit of finite-dimensional lattice inte- 
grals when the lattice spacing goes to zero. We want the regularized theory 
to preserve as many symmetries which the original continuous theory has as 
possible. Actually, some symmetries can be broken on the lattice (as rota- 
tional and Lorentz symmetries are) and be restored in the continuum limit, 
but it is not always straightforward and not always true that the symmetries 
are restored, indeed. For example, breaking gauge symmetry on the lattice 
would probably lead to a nonsensical, not gauge-invariant continuum theory. 
Gauge-invariant lattice action for a pure gauge theory was written by Wilson 
long time ago. But the task to put fermions on the lattice so that not only 

c Equation (hd[) is the leading chiral-noninvariant contribution in C c ff . Also terms of higher 
order in M, as well as terms of first order in M but involving derivatives of U, are allowed. 
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the gauge symmetry, but also chiral symmetry would be preserved turned out 
to be much more difficult. 

Let us define to this end Grassmann variables ip n , ip n in the nodes of the 
lattice for each quark flavor (color and Lorentz indices are not displayed) . As a 
first and natural guess, let us write the lattice counterpart of the Dirac action 
as follows 

ia 3 



+ ma A ^2,ip n ip n . (32) 

n 

We see that the action ([52]) reproduces the Euclidean action 

S ieim = - [ d**^^ - - m^} (33) 



in the continuum limit. Indeed, for free fermions 

,3 



ia' 
~2 



Expanding U ntfl + e = 1 — iaA^e^ + 0(a 2 ), we also restore the interaction 
term, and the last term in Eq. ( |32|) turns into the continuum mass term. The 
action ) is invariant under the gauge transformations when the U and ip 
are transformed according to 

Ipn -> ^nlpn , (34) 

where {f2 n } is a set of unitary matrices defined on the nodes of the lattice. 

Equation ( ^) is called the "naive lattice fermion action", and I have 
to say that, if the reader was convinced by the above reasoning that, in the 
continuum limit, it goes over to Eq. (^3|), he/she was naive, too. Our implicit 
assumption was that the fermion fields ip n depend on the lattice node n in a 
smooth manner, so that the finite difference ipn+e ~ fpn-e^ goes over to the 
continuum derivative. It turns out, however, that fermion field configurations 
which behave as ip n ~ ("l)™ 1 or VVi ~ ( — 1)" 2+ ™ 4 and change significantly at 
the microscopic lattice scale, are equally important. After carefully performing 
the continuum limit, these wildly oscillating modes give rise to 15 extra light 
fermion species with the same mass, the so-called doublers. 
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To understand it, consider first free massless fermions. Let 

(<9+V)n = ^bPn+ e>l ~ Ipn], {d~l(>) n = ~ [ipn ~ tpn-ej (35) 

be the forward and backward lattice derivative operators. The naive free mass- 
less Dirac operator is 

»L> = ~\lM+ d ») ■ ( 36 ) 
The eigenfunctions of S)f rcc are characterized by the Euclidean 4-momentum 

Up(n) = C p e ia ^ , 

where C p is a constant Grassmann bispinor. The eigenvalue equation 
5}f rcc u p (n) = — i\ p Up(n) implies 



sin(ap M ) 



Cp — —iXpCp (37) 



with 



A p = ±^y^ sin 2 (a PAI ). (38) 

(The operator ( |36| ) is anti-Hermitean and its eigenvalues are purely imaginary.) 
When ap p <C 1, we reproduce the continuum massless fermions with the 

spectrum X p — ±^/p^. Each eigenvalue ( |38| ) is doubly degenerate due to 

2 possible polarizations. The cigcnfunctionss with negative X p are obtained 
from the ones with positive X p by multiplication by 7 s . 

Note, however, that the lattice Dirac equation ( |37| ) has an additional 
discrete symmetry (Z2) 4 : for any eigenfunction u p , the function Q p u v = 
7^7 5 w p +(7r/a)e fl ( no summation over /i ) is also the eigenfunction of 5?f ree with 
the same eigenvalue X p . The operators Q p commute with 33f ree and anticom- 
mute with each other 

The functions 

Up, Qp,u p: Q[^Q V ]U P , Q[hQvQx] u pi Q[i_iQvQ\Q P }iip (39) 

form a degenerate 16-plet. 
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In the free case, each eigenstate of the naive Dirac operator is not just 16- 
fold, but 32-fold degenerate due to polarizations. In the interacting case (on a 
generic gauge field background), polarization is not a good quantum number, 
but the 16-fold degeneracy (|39| ) still holds. The naive lattice Dirac operator in 
Eq. (p2|) which can be written in the form 3}° = — §7^(2? + + 2?~), where 

(£>-^)„ = ^(^-Vn-e^n-eJ (40) 

are the covariant lattice forward and backward derivatives, still enjoys the 
symmetries 

Q»- ipn — ► (-l) n "7p7Vn (41) 

(no summation over /i). We see that if ip n changes smoothly from node to 
node, its 15 doublers (^) wildly oscillate on the microscopic lattice spacing 
scale. We might call these modes "unphysical" but they would not listen to 
us and contaminate with a vengeance any numerical lattice calculation we 
might wish to do. Some way to get rid of them should be suggested, otherwise 
QCD, the theory involving only 6 quarks with different masses, would not be 
operationally defined. 

The problem is that it is not so simple. Let us look for some other lattice 
Dirac operator D ^ D° satisfying the following 4 natural conditions: 

1. At distances much larger that the lattice spacing a, D — ► — Hp giving rise 
to a massless fermion in the continuum limit ^ 

2. All the modes of D not associated with the latter are of order 1/a (no 
doublers!). 

3. D is local. In other words, the matrix elements D nn ' decay exponentially 
fast at large distances \n — n'\ ^> 1. 

4. Chiral symmetries (||), (|3|) of the massless fermionic action are not broken 
by the regularization explicitly [the singlet axial symmetry (j^) is eventu- 
ally going to be broken due to noninvariance of the fermionic measure, 
but we require the absence of explicit breaking in the regularized La- 
grangian]. This seems to imply the condition 'Dj 5 + j 5 T) = 0. 

d Adding a finite mass term to £> presents no difficulties [sec Eq. . 
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The no-go theorem due to Nielsen and Ninomiya tells us, however, that 
such S does not exist. To understand it, consider first the free fermion case. 
The momentum p^ is then a good quantum number, and the Dirac operator 
in the momentum representation has the form 

2>(p) = lMp)+G(p) . (42) 

The condition 4 tells us that G(p) = 0. The condition 1 implies that F^{p) = 
+ 0(ap 2 ) for ap,j, <C 1 . Now, F M (p) is a periodic function of its four 
arguments p^ with the period 2n/a. It realizes thus a smooth map T 4 — > R 4 , 
where R A is the tangent space. A look at Fig. [j] can convince the reader that 
the point on the tangent space where it touches our torus has at least one more 
pre-image. And his/her intuition would not betray him/her: this statement 
can be proven in a rigorous mathematical manner. Basically, it follows from 
the fact that the degree of the map T d — > R d is zero which means thatB 

sign [det||&,F„(p)||] = . (43) 

pre — images 
of P 

As the Jacobian of the mapping p^ — > F^(p) is equal to 1 at the point p M = 0, 
Eq.(ff3"l) implies that some other pre-images of zero, i.e. some other solutions of 
the equation system F^(p) = should be present (one can have just one extra 
solution as in Fig. |l] or more: for the "round upright torus" F fl (p) = sin(ap M )/a, 
there are 2 d — 1 extra solutions). And that means the presence of doublers in 
contradiction with condition 2. 

The only remaining possibility is that F^(p) are not continuous. Be- 
sides being ugly, it contradicts also condition 3: the matrix elements D n n' = 
D(n — n') = 7 M F M (n — n') present actually the Fourier coefficients of the pe- 
riodic function ^^F^p). If the latter is discontinuous, the Fourier coefficients 
cannot decay faster than X/\n — n'\ (otherwise, the Fourier series would con- 
verge uniformly on the torus p^ £ [0, — ) and the sum of such a series would 
be continuous). 

We have proven that a lattice Dirac operator satisfying the above four 
conditions cannot be found for free fermions, but that also means that it can- 
not be found in QCD: any D with this property should also enjoy it for any 
smooth set of the link variables^] and, in particular, for the set t/ oac h link = 1 
corresponding to free theory. 

e We believe that, in the continuum limit a — > , the characteristic fields {U} contributing to 
the path integral can be gauge transformed into the form where [7 = 1 + Ola) for all links. 
Note that the statement that characteristic fields are always smooth is just wrong: there are 
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a) (1=1 



b) d=2 



Figure 1: Nielsen— Ninomiya theorem. P and P' are different zeros of the lattice Dirac 
operator. 

2.2 Ways to Go. The Ginsparg- Wilson Way. 

If we still want to build up a lattice version of QCD, we have to relax at least 
one of our four conditions. Conditions 1 and 2 are, however, indispcnsiblc: 
a lattice theory where they do not hold just has nothing to do with QCD. 
Therefore, either locality or chiral invariance of the lattice action should be 
abandoned. 

One of the possible procedures is that only one mode of each degenerate 
16-plet of D° is taken into account in the fermionic determinant and in the 
spectral decomposition of fermion Green's functions 



instantons which, in the singular gauge, involve singularities of the gauge potential ~ x^/x 2 
at the instanton center. Note also that a gauge transformation removing the singularity at 
x = moves it to some other point or to infinity. For the Euclidean torus that means that 
we cannot simultaneously require {7 a n n n k s = 1 + 0(a) and periodicity of U. 
To the best of our knowledge, a rigorous proof of this crucial assumption is absent, but it can 
be justified by arguing that the action of field configurations which are "essentially singular" 
(so that the singularity cannot be removed by a gauge transformation) would be infinite in 
the continuum limit. 




(44) 
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etc, where Uk{n) describe the fc-th eigenmode of D° as a function of the node. 
This amounts to choosing the lattice Dirac operator in the form (2) ) 1 / 16 , which 
is not local. A similar method is sometimes used in the practical lattice cal- 
culations, but besides purely technical inconveniences it is unsatisfactory from 
a philosophical viewpoint: we would like to have a local lattice approximation 
for a local field theory. 

But then the chiral invariance (|5|) is necessarily lost. Though renouncing 
chiral invariance is also not desirable - when regularizing the theory, we should 
try to preserve as much of its symmetries as possible - it is still considered as 
the least of evils. 

Two ways of chiral noninvariant lattice regularization have been known 
for some time and used in practical calculations: (i) Wilson fermions and (ii) 
Kogut-Susskind or staggered fermions. We will describe here the first method 
which consists in adding to D° the term ~ aD~'Dt = a2?+2?.T, where the 
covariant lattice derivatives are defined in Eq. (|4C|). Thus, the Wilson- 
Dirac operator is defined as 

V W = + V-) (45) 

where r is an arbitrary nonzero real constant. For free fermions, T>^ — > 
and 

I" 

which is the lattice laplacian. Passing to the momentum representation, we 
obtain 

®ZM = \ sin(a^) + * Y, sin 2 (^f) . (46) 

CI (X V a / 

The second term has the form of momentum-dependent mass. For small 
JV ^ l/ a ' ^ can be neglected and the continuum massless Dirac operator 
is reproduced. In contrast to 2)°, the operator ( p6| ) is not anti-Hermitean, and 
its eigenvalues are complex. What is important is that, for not small p lll the 
absolute values of the eigenvalues of 'S) w , 

are of order 1/a. The doublers disappear. At — (2,0,0,0), the eigenvalue 
© is S:; at p M =(0,s,0,s), it is ^, etc. 
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The chiral symmetry is broken, however, and it is messy. In principle, 
when the continuum limit a — ► is taken, the effects due to the breaking of 7 s 
invariance must be suppressed, but in this particular problem the continuum 
limit with restoration of chiral symmetry is rather slow to reach, and it is even 
not shown quite conclusively that it is reached at all. In particular, it is difficult 
to render pions light. In practical calculations, it is achieved by introducing a 
large bare quark mass of order g 2 (a) /a and fine-tuning it so that the effects 
due to two chiral noninvariant terms - the Wilson term and the bare quark 
term - would cancel each other. Needless to say, this is a rather artificial and 
unaesthetic procedure. 

As we see, this Nielsen-Ninomiya puzzle defies attempts to solve it. A 
recent remarkable observation is that the best strategy here is to follow the 
example of the Alexandre the Great and just cut it through! The adequate 
sword was forged back in 1982 by Ginsparg and Wilsorfl. They suggested to 
consider the lattice Dirac operators satisfying the relation 

7 5 D + £>7 5 = a®7 5 £) . (48) 

The anticommutator {D,j 5 } does not vanish which means that the lattice 
Lagrangian is not invariant with respect to chiral transformation 

5tp n = iaj 5 ip n , Srp n = iarp n ~/ 5 , (49) 

a lattice Euclidean counterpart of Eq. (Q). 

It took 16 years to realized that the lattice fermion action 

S F = a 4 ^7/i n D„„'Vv (50) 

nn' 

(color and spinor indices being suppressed), with £> satisfying the relation (Eq), 
is invariant with respect to the following transformations: 



= zcry 5 
Sip — iaip 



1 - i«s) 

1-iaS 



4, 

7 5 . (51) 



If 2) is local (in the sense of condition 3 in the Nielsen-Ninomiya list), Eq. (51) 
is as good a lattice approximation of the continuous chiral symmetry (|^) as the 
trivial (^) . In particular, the pions would automatically be light (massless in 
the chiral limit), and no fine tuning is required. But condition 4 above is no 
more satisfied and one can hope now to find a local D not involving doublers. 
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The problem is still not trivial: as we will see a bit later, many solutions of the 
Ginsparg- Wilson relation ( pl8| ) can be found and most of them are not wh 
we are looking for. The simplest good solution was suggested by Neuberger 
It has the form 



1 



.4 



(52) 



where A = 1 — a1) w and T) w is the Wilson-Dirac operator ( |45| ) with r > 1/2. 
In particular, for r — 1 and for the free fermions, we have 



aS(p) = 1 - 



1 



2 V sin 



2 ( apu 
2 



7 Al sin(ap M ) 



l + 8£_sm 2 (^)sm 2 (^) 



1/2 ' 



(53) 



The eigenvalues of ( ]53| ) arc different from zero provided p^ ^ . In particular, 
for 

-0,0,0 , p M =(-,-0,0), p M = , 

Va / \a a / V a a a I 



and 



/7T 7T 7T 7T\ 

V a ' a ' a ' a / 



the eigenvalues — iX of £> are all equal to 2/ a. The doublers are absent. A 
second look at Eq. (|5|) reveals a beautiful feature displayed in Fig. |[ the 
eigenvalues of !D lie on the circle 



(Re XY 



Im A 



(54) 



This property holds also in the interacting case. Using the property 33 1 = 
7 5 £>7 5 and the Ginsparg- Wilson relation (Bq), it is not difficult to see that the 
operator V = 1 — aD is unitary, i.e. its eigenvalues lie on the circle {e 1 ^}. 
The eigenvalues of D = (1 — V)/a he on the circle in Fig. 0. 

The function (^) is analytic on the torus € [0, — ) which means that 
its Fourier image decays exponentially at large distances. The Dirac operator 
thus constructed is local. In the interacting case, it stays local if the gauge 
field is smooth enough, i.e. the link variables U n ,n+e are sufficiently close to 
1. 

As was mentioned, the singlet axial symmetry (^) is anomalous which 
shows up in the noninvariance of the fermionic measure. The measure 
Yin dipndipn is obviously invariant, however, with respect to the ultralocal 
transformations (p9|). This follows from the fact that Tr{7 5 } = 0. On the 
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Im 




Figure 2: The circle of eigenvalues for Neuberger's operator. The eigenmodes with the 
eigenvalues marked by crosses are related by a 7 s transformation. 



other hand, Eq. (|l^) relates the modification of the measure to the operator 
trace of 7 5 , 



d A x 



k 



(55) 



which is thus zero too. For the naive Dirac operator, only the zero modes 



contribute to the sum (55) which means that the number of the right-handed 
and left-handed zero modes of D° should be equal. That also follows from 
the fact that 23° commutes with defined in Eq. (^|). It is not difficult to 
see that if u k is, say, right-handed, Q[^Q u ]Uk and Q\^Q V Q\Q P ]Uk are right- 
handed too, but Qfj.Uk and Q\nQ v Q\]Uk are left-handed. Thus, instead of a 
single right-handed zero mode in (the lattice approximation for) the instanton 
background, we have a degenerate 16-plet with 8 right-handed and 8 left- 
handed modes, which are no more necessarily zero modes. The vanishing of 
the index of S° is closely related to the identity (|43|). Indeed, the sign of the 
Jacobian det Hc^i*},^)! describes the orientation of a neighborhood of the pre- 
image Pi £ T 4 with respect to the orientation of the tangent space R A onto 
which it is mapped. This orientation is obviously related to chirality. 

The absence of the anomaly is one of the diseases of the naive lattice 
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Dirac operator. For the operators of the Ginsparg- Wilson type, the situation is 
different. Chiral symmetry is now implemented as in Eq. ( pl| ) and, generically, 
the measure is not invariant with respect to these transformations. We have, 
instead of Eq. ( |l5| ) 



Even though r Tr{7 5 } = (we derived this in the basis involving the eigenvalues 
of D°, but it is true in any basis), Tr {7 5 a£>} need not vanish and the anomaly 
is there. 

It is instructive to see how a nonzero operator trace in the right-hand side 
of Eq. (|56|) is obtained in the basis involving the eigenvalues of D. First, it 
is still true that for any eigenfunction Uk of J), j 5 Uk is also an eigenfunction 
with the eigenvalue A' fe = —\* k (complex conjugation appears and this is what 
distinguishes the Ginsparg- Wilson case from the continuum or naive lattice 
Dirac operators). Thus, for almost all eigenstates on the circle in Fig. |J), 
Uk and 7 5 Wfe have different eigenvalues, are orthogonal to each other, and the 
corresponding contribution vanishes. The only exception are two points on the 
circle: A = and A = 2i/a where A& = Al, and the eigenstates can have a 
definite chirality. But the doublers — iX = - obviously give zero contribution 
in Eq. (|5^). Qnly the zero modes of 2) are relevant. We have derived the lattice 



The right-hand side of Eq. ( p7\ ) presents a functional depending on the link 
variables {U}. By definition, it presents a sum over all lattice nodes of some 
local expression. After some work, one can be convinced that it goes over to 
the topological charge ( |l7j ) in the continuum limit. 

3 Continuum Theory: Some Exact Results 

We abandon the lattice now and will discuss in this section various aspects of 
chiral symmetry in the continuum limit. Sometimes, we will think in terms of 
quarks, of symmetries @ and (||), and of the order parameter ( |23| ) associated 
with the spontaneous breaking of the flavor-nonsinglet symmetry. Sometimes, 
we will describe the system in terms of the pseudo-Goldstone degrees of freedom 
and the effective Lagrangians (|26|) , (pi|) , and (|3(]). Sometimes we will confront 
the two languages using the philosophy of the quark-hadron duality — that is 
how most of the results discussed in this section, a bunch of beautiful exact 
theorems of QCD, will be obtained. 




(56) 



index theore 




(57) 
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3.1 QCD Inequalities. The Vafa-Witten Theorem 

As was discussed above, the octet of pseudoscalar mesons (71", K, if) can be 
interpreted as that of the pseudo-Goldstone particles appearing due to the 
spontaneous chiral symmetry breaking according to the pattern (^2|) in the 
massless limit. That is the reason why the pseudoscalar mesons are lighter than 
those with other quantum numbers. It is interesting that the latter statement 
can be formulated as an exact theorem of QCD without any reference to the 
(experimental!) fact that the chiral symmetry is broken. 

Consider a QCD-like theory with at least 2 quark flavors and assume that 
these quarks (denote them by u and d) have equal masses m u — = m. 
Consider a set of Euclidean correlators 

C r (x,y) = (r d (x)J du (y)) vac , (58) 

where J ud (x) are flavor-changing bilinear quark currents J ud — uTd with Hcr- 
mitean 

r = 1, 7 5 , «7 M , 7 M 7 5 , io^ . 
At large distances, the correlators ( p8| ) decay exponentially 

C r (x,y) oc exp{-Mr|x-2/|} (59) 

where Mr is the mass of the lowest meson state in the corresponding channel^] 
On the other hand, the correlators (|5^) of the quark currents can be pre- 
sented in the form 

C T (x,y) = -Z- 1 J d^ A Ti{Tg A (x,y)Tg A (y,x)} , (60) 

where 

d ^ A = n^(^ndet||m / -^||exp|- i L j F^F«„d 4 x} (61) 

is the standard QCD measure and Q A {x,y) is the Euclidean Green function 
of the u— and d— quarks in a given gauge-field background. Note that, when 
writing down Eq. (|6^) , we used the fact that J ud is not a singlet in flavor [oth- 
erwise, the disconnected contribution oc Tr {TQ A (x, x)} Tr {TQ A (y, y)} would 
appear in the right-hand side] . In addition the assumption m u — was made 
[otherwise, we would have two different Green's functions G A (x, y) ^ Q A {x, y)\. 

assume here that the quarks are confined, otherwise the whole discussion is pointless. 
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An important nontrivial relation 

7 5 £4(z,yb 5 = 0\(y,x) (62) 

holds. To understand it, write the spectral decomposition for g A (x,y), 

u k (x)u\(y) 



g A {x,y) = (ijj(xWj(y)) A = 



- «Ai 



(63) 



(cf. Eq. (fl4|)). Using the symmetry u k — > 7 5 itfc, Xk —> —Xk, we obtain 

k 



u p( x ) u l(y) 



m + iX v 



' m — iX„ 



t 



= ^(y,a:), (64) 



as annonced. We see that the pseudoscalar correlator 

~(Tr{7 5 g A (x,y)7 5 g A (y,x)}) = (Tr{\g A (x,y)\ 2 }) 

plays a distinguished role - it presents an absolute upper bound for any other 
such correlator. The fastest way to show this is to expand the 4x4 matrix 
g A (x,y) over the full basis 



g A {x,y) = s{x,y) +7 p(x,y) +i-y^v tll {x, y) + 7^7 a^(x,y) 
1 

— 1 

2 



(65) 



Then 



\c${x,y) = Jtt{ 7 \m (,-.,, )-\7 !('/..<■ ij- 



^{^(^y)! 2 } = | S | 2 + W 2 



1^1 + ' 2'"^ 



olVl 2 , (66) 



but, say 



-^Cfts,!/) = ^Tr{g A (x,y)g A (y,x)} = ^Ti{g A (x,y) 7 5 g\(x,y) 7 5 } 



1 



N 2 + bl 2 -KI 2 -KI 2 + ^. 



(67) 
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The inequalities 



\C$(x,y)\>\C£{x,v)\ (68) 

in any given gauge background, the positivity of the measure (^3l]) in Eq. (|60|), 
and the asymptotics imply that the mass Mps of the (lightest) pseu- 
doscalar meson in the ud channel is less or may be equal to the masses 
Ms, My , Ma, Mt of the lightest scalar, vector, axial, and tensor states.^ 
Let us emphasize again that this statement can be justified only in the 



theory with the positive measure (61) [e.g. in the theory with nonzero vacuum 
angle 0^0, the measure is not positive and pseudoscalar states need not to 
be the lightest], with equal quark masses, and only for the those states that 
are not flavor-singlet. For flavor-singlet states it need not be true. Consider, 
for example, the theory with just one quark of a large mass. Then the lowest 
meson states would be made of gluons and would know nothing about quarks. 
The lowest glueball state is believed to be scalar rather than pseudoscalar. 

In Sec. 1, we have mentioned already the Vafa-Witten theorem saying that 
the vector isotopic symmetry is not broken spontaneously in QCD. Now we are 
ready to prove it. Indeed, if such a breaking occured, the massless Goldstone 
scalar particles would appear in the spectrum. The inequality Mps < Ms 
implies that a massless pseudoscalar particle would also exist. But the theory 
with m u = ma ^ (which duly enjoys the exact isovector symmetry, a possible 
spontaneous breaking of which is under discussion now) has no exact axial 
isotopic symmetry, and there are no reasons for the massless pseudoscalar 
state to exist. So, it does not exist, hence the massless scalar does not exist 
either, and the isovector symmetry is not broken. 

Many more inequalities of this kind (e.g. Mn > M,r or M 7T + > M^o) can be 
formulated, but their proof relies on some extra assumptions. We address the 
reader to Ref. 10 for a nice recent review (see also S. Nussinov's contribution 
in this Volume). 



3.2 Euclidean Dirac Spectral Density 

Consider the Euclidean Dirac operator Tp in a given gauge field background 
A^(x). We assume that the system is placed in a finite 4- volume so that 
the spectrum oilp is discrete. Let {\k} be the background-dependent set of 

9 We were a little bit sloppy here. The inequalities dr3s| ) hold strictly speaking only for 
tree correlators, not for renormalized ones. However, renormalization only brings about 
multiplicative factors, which do not depend on distance. Thus, taking \x — y\ to infinity 
before the limit Atjv — > oo is done, we can ensure that the inequalities (pa) for renormalized 
correlators at large distances are fulfilled. 
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eigenvalues of Tp. The spectral density^ is denned as follows: 



(69) 



where the average is done with the weight function (J6l|) . The 7 s symmetry of 
the spectrum implies that p(A) is an even function of A. 

In contrast to solids or nuclei, the spectral density (|6|) is defined in the 
Euclidean space and seems to have no direct physical meaning. There are, 
however, a set of remarkable identities which relate the spectral density of the 
Euclidean Dirac operator to physical observables. The simplest such identity 
relates the spectral density at "zero virtuality" A = to the quark condensate. 

To derive it, set x = y in the spectral decomposition (|63|), integrate it over 

_!_,/■'-,■ - i 1 i , 1 1 w 11 - f / 11 -1 1 1 til;-, .irnl'-wrlllii ■ \ \ ■ / \ 1 - 1 1 I i \ fr-illnvi I "1 . . 1 , 1 ^ i i ■ 1 i 1 i ir/llirlit / ( i I \ Xll 

P(A), 



d x, and perform the averaging over the gauge fields with weight (|6 

using the symmetry p(—X) 



view of the definitions (|23|), (|2j), and (|69|) 
and assuming the reality of E, we obtain 



rn 



p(X)d\ 
m — iX 



2m 



p{X)dX 
A 2 + m 2 



(70) 



To understand better this formula, let us look first what happens for free 
fermions. As there is no physical dimensionfull scale in this case [remember 
that Afe in Eq. ( |69[ ) are eigenvalues of the massless Dirac operator], p(X) = CA 3 
on dimensional grounds. By counting the eigenvalues of the free Dirac operator 



A(tv) 



2tt 

T 



\ 



E 



(71) 



[antiperiodic boundary conditions for the fermions in all 4 directions are chosen, 
are integer, and each level ( fFH ) involves an extra 2A^ c -fold degeneracy] in 
the 4D ball 1/L < A < A, it is not difficult to determine C = A c /(4tt 2 ). Thus, 



(72) 

''The notion of spectral density and the definition (^) are also widely used in condensed 
matter physics and nuclear physics. It is especially useful if a system is disordered or involves 
elements of disorder like it is the case for electron spectra in most solids or for energy 
levels in complicated nuclei. It makes sense also for ordered systems (such as metals). 
In this case, rather than averaging over stochastic external field, one averages over some 
interval of eigenvalues AA much larger than characteristic level spacing, but much less than 
a characteristic scale of A on which p(A) is essentially changed. 
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In the interacting theory, the spectral density behaves as p(A) cx A 3 for A 
much greater than the characteristic hadron scale /ihadr so that interaction is 
weak. To be more precise, the power A 3 can be multiplied by an anomalous 
dimension factor 

A 



hi- 



/ihadr 

A recent one-loop calculation0 implies that a is nonzero and negative. 



We see that the integral in Eq. (70) diverges quadratically in the ultra- 



violet. The same result can be obtained directly by calculating the fermion 
bubble graph in the momentum representation 

<*w°» Km ^- (73) 

Thus, strictly speaking, the formula (173) does not make much sense as it stands. 



Note, however, that even though the (purely perturbative) contribution (73) 
diverges in the ultraviolet, it vanishes in the chiral limit to — > 0. The whole 
point is that in QCD, the integral (^) acquires an additional nonperturbative 
contribution coming from the region of small A which survives in the "contin- 
uum chiral thermodynamic limit" {first V — > oo, then to — > 0, and only then 
the ultraviolet cutoff is lifted Ajjv — > oo). The fact that chiral symmetry is 
broken spontaneously means that the vacuum expectation value (tp (O)"0(O)} is 
nonzero in this particular limit. 

Obviously, the necessary condition for the condensate to develop is p(0) ^ 
0. Neglecting all terms which vanish in the continuum chiral thermodynamic 
limit defined above, we obtain finally the famous Banks-Casher relatione?! 

(^(0)^(0))vac = E = 7rp(0). (74) 

Note that the result does not depend on flavor which tells us again that the 
flavor vector symmetry is not broken. 

Not onlx p(0), but also the form of p(A) at small A -C /ihadr can be 
determinedEj Consider the theory with Nf > 2 light quarks of common mass 
to. Let us study an integrated correlator in this theory, 

d 4 x(S a (x)S b (0)) = 1 J d i xd 4 y(S a (x)S b (y)) , (75) 

where S a (x) = ■ijj(x)t a ip(x) and t a is the generator of the SU(Nf) flavor group. 
Fix a particular gluon background and define 

C ab \ A = ~ J d A xd A yTi {t a G A (x,y)t b g A (y,x)} . (76) 
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Figure 3: Pseudogoldstone loop in scalar correlator. 



Substitute here the spectral decomposition ( p3| ) for Ga(x, y), do the integration 
and perform averaging over the gluon fields trading the sum over eigenvalues 
for the integral over the spectral density (|69|). We obtain 

S ab f°° P (X)dX 




2 J_ x {m-i\y 

- S°»r ^-*P d X, (77) 
J Q (m 2 + A 2 ) 2 

where the property p{— A) = p(A) was used. 

On the other hand, the same correlator can be saturated by physical states 
among which the pseudo-Goldstone states play a distinguished role. Consider 
the 1-loop graph in Fig. || describing the contribution of the 2-Goldstone in- 
termediate state ~ J (0\S a \4> c (/) d }((l) c (l) d \S b \0} (obviously, the one-particle state 
does not contribute because the pseudo-Goldston mesons are pseudoscalars 
while S a (x) is scalar). To calculate it, we need to know the vertex (0\S a \(f> c (j) d ) 
which can be determined via the generating functional of QCD involving scalar 
sources u a coupled to the current S a . Adding the source term u a S a to the La- 
grangian amounts to adding u a t a to the quark mass matrix M.. The latter 
enters also the mass term ( |30| ) in the effective Lagrangian. Expanding U up 
to the second order in <p a and varying it with respect to u a , we obtain 

(O|5 Q |0 c d ) = -^d abc . (78) 
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The vertex is nonzero only for three or more flavors. 

Now we can calculate the graph in Fig. Actually, we cannot because the 
integral diverges logarithmically in the ultraviolet, but anyway the effective 
theory is not valid at high momenta (technically, the divergence is absorbed 
into local counterterms of higher order in p char and m). Only the infrared- 
sensitive part of the integral is relevant. A simple calculation gives 



32^ N f \F*J MLlr 



Compare it now with Eq. (77). Note first of all that the constant part p(0) 
does not contribute here, 

\2 

' dX = . 



(m 2 + A 2 ) 2 

Thus, only the difference p(A) — p(0) is relevant. It is easy to see that, in order 
to reproduce the singularity ~ In Mj ~ lnm, we should have p(A) — p(0) — C\X\ 
at small |A|. Substituting it in Eq. (|77|) and comparing the coefficient of lnm 
with the coefficient of In in Eq. ( |79[ ) , we finally obtainO 

PW = ~+ TzZs^lps) W +o(X). (80) 
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Thus, for Nf > 3, the spectral density has a nonanalytic "dip" at A = 0. The 
behavior is smooth in the theory with two light flavors. Physically, it is rather 
natural that the more is number of flavors, the stronger is the suppression of 
p(0). The determinant factor in the measure (|6ll ) punishes small eigenvalues, 
and the larger is Nf, the more important is this factor. By "analytic continu- 
ation" of this argument, one should expect a nonanalytis bump rather than a 
nonanalytic dip at A = in the case Nf = 1. Indeed, Eq. (^) displays such a 
bump. One should not forget, of course, that the whole derivation was based 
on the effective chiral Lagrangian approach and does not directly apply to the 
case Nf — 1. Some additional, more elaborate reasoning shows, however, that 



a bump at Nf = 1 as predicted by Eq. (80) is thereBJ-.lt is also confirmed by 



a numerical calculation in the instanton liquid modelta 
3.3 Chiral Symmetry Breaking and Confinement 

Look again at Eq. (||). The axial current entering the left-hand side is an 
external current in a sense that no dynamical field is coupled directly to j^. 
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But the fields entering on the right-hand side are dynamical gluon fields present 
in the QCD Lagrangian. 

In the chiral (left-right asymmetric) gauge theories like the standard elec- 
troweak model, both vector and axial currents are coupled directly to the phys- 
ical gauge fields. Anomalous divergence of such current would mean explicit 
breaking of the gauge invariance which is not nice. Therefore, in chiral the- 
ories, one should always take care that such purely internal anomalies would 
cancel out at the end of the day. In the Standard Model, they do. 

Let us discuss, however, purely external anomalies in QCD which are not 
related to breaking of any symmetry but just mean that certain correlators 
involving external currents are singular .[] 

As a simplest nontrivial example, consider the theory with two massless 
flavors and look at the correlator 

H (q) = * J TOtoj?(Q)}>7< e^d^x , (81) 

where A, B are flavor indices and Ti is the external flavor-singlet "magnetic 
field." Fl The correlator (|Tj) is nothing but a three-point vacuum expectation 
value (ufi in kinematics in which one of the external momenta associated with 
the vector current is set to zero. 

The one-loop calculation of the corresponding graph displays a singularity, 

H {i) = -2L q -^- ■ N c ■ h*» . (82) 

p 2-k z q A 2 

The imaginary part of this amplitude.-is also singular, ~ 5{p 2 ), which can be 
related to the masslessness of quarksll3 However, the quarks (in contrast to 
electrons in QED) do not exist as physical particles due to confinement and 
one can ask where does the singularity in the imaginary part Im Kj~' n (q) oc 
5(q 2 ) come from? This is a good question, the answer is better still: the 
singularity ~ 1 j q 2 comes from the propagator of a massless Goldstone boson, 
which appears due to the spontaneous chiral symmetry breaking and which is 
directly coupled to the axial current (see the middle graph in Fig. ||). 

Let us ask now: can one reproduce the singularity in Eq. ( p2| ) without 
Goldstone bosons and without spontaneous chiral symmetry breaking, but in 
some other way ? 

8 The conventional anomaly (^) is related to the correlator involving both internal (j™) 
and external currents and can be called mixed in this setting. 

J The quotation marks distinguish "H which couples to the baryon charge from the physical 
magnetic field which has the matrix structure diag(2/3, —1/3) and is a mixture of isotriplet 
and isosinglet. But we are not interested in dynamics of electromagnetic or weak currents 
here. In QCD proper, all color-singlet currents are external. H is just a source of such vector 
flavo-singlet current. 
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Figure 4: Saturating the external anomaly. 



As far as the theory with two light quarks is concerned, the answer is pos- 
itive: the singularity of the correlator above can, be reproduced, in principle, 
if massless baryons are present. Proton and neutron represent, like quarks, a 
flavor SU(2) doublet. There are N c = 3 quark doublets and only one baryon 
doublet |P) = \uud) and \N) = \udd). The absence of the overall N c factor 
is compensated, however, by the fact that the baryon charge of nucleon is 3 
times larger than that of the quark, and the vertex involving the "magnetic 
field " TL is 3 times larger for baryons. 

Thus, this purely algebraic anomaly matching argument due to 't HooftS 
does not rule out a dynamical scenario where the physical spectrum in the 
theory with just two massless quark flavors would not involve massless pions, 
but, instead, the massless proton and neutron. It is rather remarkable that, 
in the theories with Nt > 3, the scenario with massless baryons is ruled 
out. Suppose that, instead of the octet of massless Goldstone fields, we have 
an octet of massless baryons. The contribution of the corresponding triangle 
graph in would have the same structure as in Eq. (^) , but with the factor 

Tr {T a T b } = C 8 S ab (83) 



instead of 5 ab /2, where T a are the flavor generators in the octet representation. 
To find the Dynkin index of the octet representation Cs, it is sufficient to 
assume that a = b = 1 and decompose the octet with respect to the SU(2) 
flavor subgroup: 8 = 3 + 2 + 2 + 1. The contribution of each doublet to C% 
is 1/2 and the contribution of the triplet is 2. Adding it together, we obtain 
C*8 = 3 7^ 1/2 and the required result ( |S2"|) is not reproduced. Also a massless 
decuplet and all other possible color-singlet baryon representations would give 
the coefficient in front of the singularity much larger than that in Eq. (|si| ) , and 
the anomaly matching condition would not be fulfilled. Therefore, massless 
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baryons do not exist .[] 

We have arrived at a remarkable result. In QCD with three massless 
quarks, the assumption of confinement allowing the existence of only colorless 
states in the physical spectrum and the anomaly matching condition lead to the 
conclusion that massless Goldstone states must appear and chiral symmetry 
should be spontaneously broken. If the latter is not true, the only possibility to 
saturate the anomaly is to assume that massless quarks still exist as physical 
states in the spectrum and there is no confinement! 

In the real world, confinement and spontaneous chiral symmetry breaking 
in the limit of massless u, d, and s quarks arc experimental facts. Whether 
or not these phenomena take place in hypothetic theories with Nf > 4 is an 

open question. It is quite possible that starting, say, from Nf = 6, the small 
eigenvalues in the Dirac operator spectrum are punished by the determinant 
factor so strongly that p(0) = and, in view of the Banks-Casher relation <y4\), 
the quark condensate vanishes, and the symmetry is not broken. The anomaly 
matching argument tells then that there is no confinement in this case.| 

We have called this result - that the chiral symmetry breaking and confine- 
ment go together - remarkable. It is also somewhat misterious. Even though 
we do not understand well dynamical reasons for confinement to occur, we 
still expect that the same mechanism which works for the theory with Nf — 3 
works also for the theory with Nf — 2. But 't Hooft's argument works only 
for N f > 3... 



4 Mesoscopic QCD 

To make the path integral finite-dimensional and the spectrum of the Dirac 
operator discrete, not only an ultraviolet lattice regularization should be intro- 
duced, but one also has to consider the theory in a finite volume. If the volume 
is very large, much larger that the pion Compton wavelength, the dynamics 
of the theory is basically the same as that in the continuum limit. However, 
in practical calculations characteristic lattice sizes are not too large, and finite 
volume effects are essential. It is therefore important to estimate theoretically 

k To be quite precise, one could, in principle, saturate the anomaly with several baryon 
multiplets with positive and negative baryon charges. This possibility is so unaesthetic, 
however, that it can be rejected by that reason. 

; We know, of course, that if the number of the quark flavors is very high Nr > 16, the 
asymptotic freedom is lost and we cannot expect confinement. We are almost sure that 
quarks and gluons are not confined at Nf = 16 or Nr = 15, in which case the theory is 
asymptotically free, but has-an infrared fixed point at a small value of ct s , and the coupling 
constant never grows largely There are some reasons to believe, however, that confinement 
is lost at a smaller value of Nf, not just Ne = 15. 
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these effects. In this sectkp, we derive a set of exact relations for QCD placed 
in a finite Euclidean boxE3 They rely heavily on the presence of the chiral sym- 
metry and on the notion of quark-hadron duality used in the previous section. 
We will always assume that the size of the box L is large L ^> /ij^dr' but the 
requirement L 3> M ~ 1 is not imposed. 



4-1 Partition Function: Nf = 1 

Consider first the theory with just one light quark flavor. We assume that the 
quark mass m is complex so that the mass term in the Lagrangian has the 
form 

Cm -1 = mipRipL + rhip L ip R , (84) 

and the vacuum angle 9 is nonzero. The partition function of the theory can 
be written in the form 



Z(Q) = exp{—Ve(9, m, L)} 



(85) 



where V — and e(9, m, L) tends to the vacuum energy density e vac {9, m) in 
the thermodynamic limit L — > oo. It is important that the finite volume cor- 
rections e(#, m, L) — e vac (#, m) are exponentially suppressed cx exp{— ^hadri}- 
Indeed, the corrections appear due to modification of the spectrum of excita- 
tions in the finite box by the Casimir mechanism. At the technical level, they 
are due to modification of the corresponding Green functions by finite volume 
effects. But in the theory with Nf = 1, all physical (meson and baryon) exci- 
tations are massive and their Green functions decay as exp{— /Zhadr-R} at large 
distances. Also the modification of the Green functions (at coinciding points) 
due to boundary effects has the order exp{— /Lthadr^}- 

The vacuum energy depends on the quark mass m and vacuum angle 9. It 
turns out that e vac is actually a function of just one complex parameter which 
is me' 9 . To see that, express Z(9) as a series 



Z{9) 



E < 

q— — oo 



n<L4«(r,x)det 



-lip + m h m 



x exp 



^ Jo ^ J dx[(Fp2] } S Pj~ iq6Zq ' 



(86) 



where Z q is given by the functional integral over the quark and gluon fields in 
the sector of definite topological charge q. Consider some particular Z q and 
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suppose for definiteness that q > 0. One of the factors entering the integrand 
of Z q is the quark determinant, 



det 



— UP + m h m 

r 2 2 



<?]Jdct 



m 

k 



m —i^k 
-i\ k m 



*Y[{\l+mm), (87) 



m 

fc 



where the product J} runs over nonzero Afc's, and the factor to 9 reflects the 
presence of q right-handed fermion zero modes as the index theorem (|l^) re- 
quires [cf. Eq. (||) written earlier for the case of real masses]. If q < 0, the 
factor fh q should be substituted by m~ q . Averaging ([87| ) over the gluon field 
configurations with a given q and substituting it in Eq. (lS6l), we obtain that 
each term in the sum depends, indeed, only on the combinations z — me 
and z (being real, the partition function and e vac cx hiZ cannot, of course, be 
holomorphic in z). 

Next we assume that \m\ -C ^thadr, that the dependence of e vac on mass is 
analytic, and expand 



e vac = eo - SRe (me* 9 ) + o(m) , 

where S is a real constant which, for 8 — and real m, coincides with the 
variation of the vacuum energy with respect to mass, i.e. with the quark 
condensate^ (s/^vac- Substituting (|8^ ) in Eq. j8^ ) and comparing it with the 
Fourier expansion (Bq), we can derive 



1 

Z q (X — dde iqe e VXm cos 8+o(m) = J^^V) (89) 

[ig(a;) = I- q {x) stands for the Bessel functions which grow exponentially at 
large real x]. 

It is very instructive now to consider two limits: (i) mT,V 3> 1 and (ii) 
mT,V ^ 1 (keeping the assumption L/Jhadr 3> 1). The limit (i) is a genuine 
thermodynamic limit. The sum ( |86| ) is saturated by many terms with large 
values of Actually, one easily derives (for any mHV) 



\<1 



2 ) vac = mXV. (90) 



m A finite nonzero value of the quark condensate was exactly our assumption when we wrote 
the expansion (p8h. 
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The quantity \ — (q 2 }vac/V can also be expressed via the correlator of topo- 
logical charge densities, 



X 



1 



167P 



;Tr{^F Mi/ }(x) 



1 



167T' 



:Tr{^F Ml/ }(0) 



(91) 



and is called the topological susceptibility of the vacuum. Taking a double 
variation of Eq.(|8^) with respect to 9, we relate it to the curvature of the 
function e vac (8) at 9 = 0: 



= — lim 



1 



d 2 Z{9) 



y^oo VZ(0) d9 2 



d 2 e 



9=0 



d6 2 



8=0 



The region 



A*hadr 



< L < (mS) 



-1/4 



(92) 



(93) 



can be called mesoscopic. The word is borrowed from the condensed matter 
physics where mesoscopic system is defined as a system involving a large enough 
number of particles for the statistical description to be adequate, but which 
is not sufficiently large for the boundary effects to be irrelevant. The same 
is true in our case. In particular, when mT,V <C 1, the partition function is 
well approximated by just the first term q = and, if we are interested in the 
mass dependence, — by the terms q = ±1 corresponding to the instanton and 
anti-instanton sectors. 

If mTiV <C 1, the fermionic condensate is provided by instanton-like config- 
urations supporting a single fermion zero mode. If mSF ^> 1, the condensate 
appears due to finite average density of small but nonzero eigenvalues of the 
Dirac operator, according to the Banks-Casher formula ([74]). Quite remark- 
ably, the value of the condensate is the same, however, and does not depend 
on the parameter mT,V. 

Note that, if m is real and 9 — (more generally, if the combination 6 p hy S — 
9 + arg(m) on which everything depends vanishes), all terms in the partition 
function ( p6[) are positive and the path integral for Z{9) has a probabilistic 
interpretation. And if not - then not. 



4-2 Partition Function: Nf > 2 

The spectrum of the theory with several light flavors involves light pseudo- 
Goldstone states. If the length of our box is comparable with the Comp- 
ton wavelength of pseudo-Goldstone particles, finite volume corrections in 
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e(9, m, L) in Eq. ( pq ) are not exponentially small and should be taken into 
account. We will be interested here in the region 



The scale ~ m" 1 ^^^ entering Eq. (|93|) and lying in the middle of the 
newly defined mesoscopic interval (^) is relevant in the multiflavor case too. 
The formulae we are going to derive in this section will be universally valid, 
however, for any value of mT,V as long as the condition (Q) is satisfied. 

The main idea is to present the partition function as the path integral in 
the effective theory (26), ( |30| ) and notice that, when the condition LM, <C 1 
is fulfilled, only the zero Fourier harmonics of the pseudo-Goldstone fields are 
relevant.^ Thus, we can disregard the kinetic term, and the path integral is 
reduced to an ordinary one. We want to calculate the partition function at 
arbitrary 9 and need to know the form of the effective potential at arbitrary 9. 
In the case Nf = 1 all physical quantities depended on the combination me , 
but not on m and 9 separately. By the same token, for Nf > 2, they depend 
on the combination Aie lS ^ Nf (where M. is a complex mass matrix). The proof 
is quite analogous and is based on the multiflavor version of Eq. (|S7j) , 



M v ^m/ihadr 



(94) 



det 



-Vp + M 



f-7 5 



M 



1 



= [det ||7W t ||] t? 



11 dct 



-iXh 



-i\k 
M 



= [dct |LM f ||] q Y[ det \\Xl + M^M\\ . (95) 



Similar to the case Nf — 1, we see that the physical vacuum angle on which 
all physical quantities depend is 



6> p hys = + arg(det||M||) 



(96) 



rather than just 9. In particular, the statement made earlier that the vacuum 
angle in QCD is either zero or very small refers to the combination (|96|). 
Thus, we can write 



Z Nf (9) = J d[/exp{vERe Tv{Me l0/N 'U^} } , 

pseudo-' 
V ^ 



(97) 



"One can observe it by looking at the pseudo-Goldstone Green function in finite volume, 

\ 



Gtt{x) 



L ) 
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where the integral is done over the group SU(Nf) with the proper Haar mea- 
sure. Integrating it further over d9 as in Eq. (|S9|), we find the partition function 
in the sector with a definite topological charge Z q . In the simple case M. — ml, 
a beautiful analytic expression for Z q can be derived, 



Z q = det 



Iq(K) Iq+l{n) 
Iq-l(K) Ig(K) 

Iq-Nj + l{^) 



lq+N f -l{n) 



(98) 



where k = \mT,V\. 

Let us assume 8 = , m real, and let us study the quantity 



(W>) 



m,V 



_L d_ 

NfVdm 



InZ , 



(99) 



the quark condensate at finite mass and finite volume. In the limit k — > oo, the 
integral in Eq. (|97|) has the main support in the region U ~ 1, and we obtain R 



e N fK 



(100) 



so that the condensate (Q) tends to a constant S as it should. 

In the opposite limit k«1, -Zq^o are suppressed compared to Zq as n^ Nf 
reflecting the presence of Nf\q\ fermion zero modes in the sector with the 
topological charge q. The condensate @) is also suppressed oc n. This is due 
to the fact that, in the multiflavor case, the condensate plays the role of the 
order parameter signaling the spontaneous chiral symmetry breaking. Strictly 
speaking, spontaneous symmetry breaking never occurs in quantum systems 
at finite volume in the sense that the true ground state wave function is always 
symmetric even though the gap separating it from excited asymmetric states 
may be very (exponentially) small.[] So , for Nf > 1, the condensate should 
vanish in the chiral limit m — > no matter how large the fixed volume V is. 
And it does. 



°Note the presence of the power-like preexponential factor in Eq. (100). It reflects the fact 
that, while the condition LM^ < 1 is satisfied and higher Fourier harmonics of the pion 
fields are decoupled, we are not in the true thermodynamic limit yet. In the latter, the 
partition function is given by a simple exponential as in Eq. (p5|), and the finite- volume 
corrections to e(8,m, L) are suppressed as exp{— M-^L}. 

p Never say never: the statement above is not correct for supersymmetric systems. But QCD 
is not supersymmetric and we may forget this subtlety. 
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4-3 Spectral Sum Rules 



Let us return to a simpler case, Nf — 
positive, and rewrite the result (P9|) as 



1. Assume that m is real and q is 




CI q ( K ) 



(101) 



where the averaging is done over the gluon field configurations with a given 



topological charge q. We can expand now the two sides of Eq. (101) in m 



and compare the coefficients of the expansion. The first nontrivial relation is 
obtained when comparing the terms cx m q and cx m q+2 . We obtain 




4(« + 1) 



(102) 



where the sum runs over the nonzero eigenvalues and the symbol ((•■■)) means 
averaging with weight exp{— S g } Y[' k A| (the second factor is the fermionic de- 
terminant in the lim it m — > 0, with the discarded zero modes). 

Expanding Eq. (101) up to the order ~ m g+4 , we obtain the following sum 
rule: 



i 



\ 2 \ 2 



16(g + !)(? + 2) 



(103) 



Similar relations can be derived in the multiflavor case. Expanding Eq. (|9^) in 
mass and comparing the te rms ~ m N f q , ~ m N f q+2 , and ~ m N f q+A , we obtain 
the same relations (102), (103), with q being substituted by q + Nf — 1. We 
can extract more information by considering the expression for the partition 
function with different quark masses. Two different sum rules can be derived 
by considering the terms ~ m Nfq+A of the expansion, 



(T,V) 4 




16 [(q + Ntf-l] 



16(q + N f )[(q + N f f 



(104) 



Subtracting one sum rule in Eq. (104) from the other and setting Nf = 1, 
we reproduce Eq. (103), but the relations (104) are valid only for Nf > 2. 
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Actually, the average 

((EfcUMfe)))? diverges at N f = 1. For N f = 2, there 
are two sum rules at each level. For Nf = 3, three different sum rules can be 
derived starting from the terms ~ rn N f q+6 of the expansion, etc. 
The sum rules ( |102| ) can be presented as the integrals 

vj 

where p(X) is the microscopic spectral density fl69|) in the sector with a given 



topological charge q.. The sum rules ( |103| ) and the first sum rule in Eq. (104) 
are not expressed in terms of p(A) only, but also via a certain integral involving 
the correlation function p(Ai, A2). The sum rule for ((Yl'k^i=ip(^/^k^f^p)))q 1S 
expressed in terms of an integral of the correlator p{Xi, X2, X3) , etc. Using 
ingenious arguments which are beyond the scope of our discussion here, J. 
Verbaarschot and collaborators managed to determine the functional form of 
p(X) and of all relevant correlators. For example^ 

p*f«W = [j 2 Nf+lql &vx) 

- J Nf+M+1 (ZVX)J Nf+lql _ 1 (VVX)] . (105) 

The expression is valid for A <C /ihadr, and the function is "alive" in the region 
of very small eigenvalues A ~ 1/SV. In the thermodynamic limit XY>V 00, 
the function (105) tends to the constant S/7r in agreement with the theorem 

Spectral sum rules are well adapted to be confronted with the numerical 
lattice calculations. The main interest here is not so much to "confirm" these 
exact theoretical results by computer, but, rather, to test lattice methods. 
This was a challenge for lattice people for some time. By now reasonably 
good numerical data were obtained, they agree well with the theorycil One can 
expect that the accuracy will grow substantially if the algorithms based on the 
lattice fermion action with exact chiral symmetry (as was discussed in detail 
in Sec. 2) are implemented. 



5 QCD at finite 9 

We know that in the real world the gauge group is SU(3) and we have 6 
fundamental quarks with some particular mass values. We also know that 
#phys defined in (^6|) is very close to zero. It presents a considerable theoretical 
interest, however, to study what happens in the same theory, but with other 
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values of the parameters. Though these "fairy" [*] variants of the theory do not 
have a direct relation to reality, in order to understand well the physics of our 
world, it is important to understand also how it varies if the rules of the game 
are changed. 

One of the ways to modify the theory is to assume a nonzero p hys ■ We 
will not assume that the quark masses coincide with their experimental values, 
but will keep them (or rather some of them) small enough for the chirality 
considerations to be relevant. The expressions for the partition function in 
finite volume as a function of 9 were written in Sec. 3. In this section, wc 
will concentrate, however, not on the mesoscopic regime m'EV ~ 1, but on 
dynamics of the theory in the thermodynamic limit mEV 3> 1. ._. 

The vacuum energy at finite 9 (its mass-dependent part) is obtainedtB by 
minimizing the effective chiral potential over U, 

e vac {0) = -E min Rc Tr {Me' 9 ^^} . (106) 

In the two-flavor case, a nice analytic expression for e vac (9) can be written for 
a generic mass matrix M, 



e vac (0) = -VyjTr{MiM} + e ig det\\M\\ + e~ je det ||7Wt|| . (107) 

In the general case it is a continuous smooth function of 9. It is, of course, 
also periodic: t vac (9 + 2ir) — e vac (9). The case of d egen erate masses is more 
subtle. Substituting A4 = ml (with real m) in Eq. ( |l07|) we obtain 



z (9) = -2mE 



cos - 
2 



(108) 



This function is still periodic in 9, but is obviously singular at 9 = tt. A 



mathematical reason for this is rather simple. The function (107) has two 
stationary points: a minimum and a maximum. At 9 = tt, they fuse together. 
The maximum and minimum of energy are plotted in Fig. g, We see that 
two smooth curves oc ±cos((9/2) cross at the point 9 = tt and, after passing 
this point, the maximum and minimum are interchanged. While e vac (9) is 
determined by the curve — 2mE cos(6*/2) at 9 < tt, it is determined by the 
curve 2m£cos(#/2) for 9 > tt. 

Physically, this is the situation of the second order phase transition: the 
energy is a continuous function of the parameter 9, but its derivative is not. 



■JThis word was coined by chess problemists. A fairy chess problem is a problem referring 
to a game with modified rules or on a nonstandard board. 
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Figure 5: Nj = 2: Minima (solid line) and maxima (dotted line) of the effective potential 
for different 9. 



A more detailed analysis shows, however, that the conclusion of the existence 
of the second order phase transition is an artifact of the approximation used. 
Indeed, for Nj = 2, the potential in Eq. (106) just vanishes identically at the 



point 9 = tt. That means that the pion excitations become exactly massless 
and chiral symmetry is restored^] 

It is clear from physical considerations that, no matter what the value of 9 
is, the symmetry ( fjlj ) is broken explicitly if the quarks are massive. For 9 = tt 
the leading term vanishes, but the breaking is still there beingjiiduced by the 
terms of higher order in mass in the chiral effective potentials The relevant 
term has the structure oc m 2 sin 2 (#/2)Tr{£/ 2 } . 

One can show that taking this into account modifies the curves in Fig. || 
in the vicinity of 9 = tt so that, at the point 9 = tt, one has two exactly 
degenerate vacuum states separated by a low but nonvanishing barrier. If 9 is 
slightly less than tt, the degeneracy is lifted, and we have a metastable vacuum 
state on top of the true vacuum. When 9 slightly exceeds tt, the roles of these 



""Another way to understand it is to look at the Gell-Mann-Oakes-Renner relation (B 
Equation ( pl| ) was originally written for 9 = 0, but a generalization for arbitrary 9 is rather 
straightforward. In our case, we even do not need it: the point 9 = n, m u = is equivalent 
to the point 9 = 0, m u = —rrid by virtue of (Bq), and we see that the pion mass vanishes. 



: J »8 



two states is reversed: the one which was stable before becomes metastablc, 
and the other way around. And this is exactly the physical picture of the first 
order phase transition with superheated water, supercooled vapor and all other 
familiar paraphernalia. 

Consider now the theory with three degenerate light flavors. By a conju- 
gation U — > RUw , any unitary matrix U can be rendered diagonal, 



U = diag(e l 



)C tf e -i(a+/J)) 



When M. = ml, the conjugation does not change the potential in Eq. ([TO 1 
For diagonal U, the latter acquires the form 



V(a,f3) = -roS 



cos a 

3 



cos 8 

1 3 



cos [ a + 8 + — 



(109) 



The function U has six stationary points, 

I: a = /3 = 0, II: a = 8 = - — , III: a = 3 = — 

3 3 



IV : a = 3 = 



20 



it, TV a : a = ~a — /?=—- 



20 



IT, 



TVb: = -a- 3 = 



20 



+ 7T. 



(110) 



The points IVa and IVb are obtained from IV by the Weyl permutations 
a <-> 8, a «-> —a — 3, etc, and their physical properties are the same. Actually, 
we have here not 3 distinct stationary points, but, rather, a four-dimensional 
manifold SU(3)/[SU(2) ® U(l)] of the physically equivalent stationary points 
related to each other by conjugation. The values of the potential at the sta- 
tionary points are 



€1 = — ottlIj cos — , en = -omLcos , 

3 3 

em = — 3mH cos — - — , eiv = mY, cos 9 . 



(Ill) 



Studying the expressions (111), and the matrix of the second derivatives of the 
potential (10E) at — tt, one can readily see that i) the points I and III arc 
degenerate minima separated by a barrier; ii) the point II is a maximum, and 
Hi) the points IV are saddle points. The physical picture we arrive at is that 
of the first order phase transition. The situation is even simpler than in the 
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Figure 6: Nf = 3: Stationary points of V(a, ji) for different 8. The solid lines are minima, 
the dotted lines are maxima, and the dashed line are saddle points 



two-flavor case since the barrier between the minima appears in the leading 
order and the subleading m 2 terms in the potential is of no concern to us in 
the case at hand. 

Figure || illustrates how the stationary points of the potential move when 
the vacuum angle is changed. At 9 = w/2, 9 = 3ir/2, etc, the metastable 
minima coalesce with the saddle points and disappear. No trace of them is left 
at 9 = 0. One can show that also for the physical values of the quark masses 
the metastable vacua are absent at 9 = 0. 

If the quark masses are not degenerate but their values are close, nice 
curves in Fig. |] are distorted a little bit, but physics remains largely the same: 
first order phase transition is robust and cannot be destroyed by small vari- 
ations of parameters. If the masses are essentially different, the phase tran- 
sition may disappear. This is true, in particular, for physical values of the 
quark masses. In this case, the vacuum energy depends on 9 analytically, as 
in Eq. (107). Unfortunately, no such simple formula for e vac (9) can be written 
for a generic mass matrix M. for three or more flavors. 

What we can do easily and quite generally is to determine the behavior of 
the function e va , c (9) at small 9 and the topological susceptibility (f)2|). Let us do 
it for arbitrary Nf still assuming, for simplicity, that M. is real and diagonal. 
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Then U can be taken to be diagonal too, 
U = diag(e Ml , ... e laN i 



N f 



For 9 <C 1, we anticipate that the vacuum values of otj will also be small. We 
can write the effective potential as 



N f 



V cS {U) = -e£] 



nij cos a, — 



j'=i 



A'/ 



= const 



AT, 



(112) 



Adding here a Lagrange multiplier term A Yl^di a j ana - minimizing the expres 



sion thus obtained over a, and A, we find 



N f ^E.iWJ 



(113) 



Substituting this result for aj in Eq. (|112|), we finally obtai 



X = m cff = gg 2 lfl=0 = yil )N S 



1 



TOjv f 



(114) 



Roughly speaking, x is proportional to the lighest quark mass (lighest quark 
masses). In actual QCD 



Xqcd 



m u m d , 



(115) 



the strange and other quarks are too heavy to be relevant. If Ai = ml, 
Eq. (114) is reducedo to x = friE/Nf. The latter result can be easily repro- 
duced by substituting the asymptotics (|10C|) for Z q in the definition 



J2 q z q 
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